Role of exchange interaction in self-consistent calculations of endohedral fullerenes 
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Abstract 

k£)Results of the self-consistent calculation of electronic structure of endohedral fullerene Ar@C6o within the Hartree-Fock and the 
7"} local density approximations are presented. Hartree-Fock approximation is used for the self-consistent description for the first time, 
"^-n It is shown that the accurate account of the exchange interaction between all electrons of the compound leads to the significant 
£T) modification of the atomic valent shell which causes the noticeable charge redistribution inside the endohedral compound. 
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1. Introduction 

Eh ■ Endohedral fullerenes represent a class of compounds con- 
C3 sisted of atoms or small molecules confined inside the fullerene 
cage. Due to the confinement of the encaged atom endo- 
O fullerenes have unique properties and a wide range of potential 
' S applications. For instance, they are used in drug delivery appli- 
cations JJ], doping of endohedral fullerenes allows one to im- 
prove power conversion efficiency of the photovoltaic devices 
l_?112|]. Quantum dot devices based on endohedral fullerenes en- 
capsulated in single-walled carbon nanotubes are reported to be 
,— ' treated ^ and proposed for quantum-information processing 

I ' Besides the potential applications the investigation of fun- 
damental properties of endohedral fullerenes is also of a great 
current interest. Since the first experimental observation of en- 
dohedral fullerenes a large number of theoretical investiga- 
tions have been performed on the photoionization process of 
the compounds within different frameworks (see, e.g., review 
|Ql and references therein). On the contrary, the experimen- 
tal studies are still not numerous what is caused by difficul- 
ties in producing sufficient amounts of purified compounds for 
the gas phase experiments. Up to now series of experiments 
have been performed on the neutral and charged endometallo- 
fullerenes O, H lU an d on noble gas endohedral fullerenes, 



OO 

o 



- 1—^ 

x 



Ar@C 60 HOI HID and Xe@C+ [LJ as well. 

In this paper we present the results of the first self-consistent 
description of electronic structure of endohedral fullerene 
Ar@C6o within the Hartree-Fock (HF) approximation. From 
atomic physics it is known [13] that an accurate description 
of the photoionization process of many-electron systems is ob- 
tained within the Random Phase Approximation with Exchange 



(RPAE) which is based on application of the HF approximation 
for calculation of the ground and excited states of the system. 
It was proven by case studies of a variety of isolated atoms and 
their ions Ill4ll and of metallic clusters as well IU5I1 . 

Considering Ar@C6o within the local density approximation 
(LDA) in [ 16] the hypothesis of hybridization of the valent 3p 
shell of the confined Ar atom was discussed for the first time. 
In the present paper we describe the hybridization phenomena 
in another way and find out the correlated effect caused by the 
hybridization. Application of the HF approximation allows one 
to take into account the exchange interaction between all con- 
sidered electrons of the system more accurately in comparison 
with the LDA. 

2. Model description 

In a number of papers it was supposed that in noble gas en- 
dohedral fullerenes A@Ceo a confined atom A is located in the 
center of the fullerene molecule Il7l Ilia, Il9ll vibrating due to 
finite temperature 12011 . In the present paper we assume that 
in case of Ar@C6o these vibrations are insignificant at room 
temperature and consider Ar atom located in the center of the 
fullerene cage. Thereby all the electrons of the endohedral sys- 
tem are moving in a spherically symmetric central field and 
one can construct the electronic configuration described by the 
unique set of quantum numbers {n, /} where n and / are the prin- 
cipal and the orbital quantum numbers correspondingly. 

In the present paper the following scheme of a shell filling 
is proposed. Electronic configuration of an isolated Ar atom is 
ls 2 2s 2 2p 6 3s 2 3p 6 . Using the atomic classification, the config- 
uration of 240 delocalized electrons in the pristine C6o can be 
written as: 
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ls 2 2p 6 3d 10 4f 14 5g 18 6h 22 7i 26 8k 30 91 34 1 0m 18 
2s 2 3p 6 4d 10 5f 14 6g 18 7h 10 , 



(1) 
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where the first 10 shells (Is ... 10m) are nodeless and the last 
6 shells (2s ... 7h) have a single radial node. Since it is com- 
monly acknowledged that fullerenes are formed from the frag- 
ments of the planar graphite sheets, it is natural to match cr- 
and 7r-orbitals of the graphite to the nodeless and the single- 
node wave functions of a fullerene, respectively 12111 . Differ- 
ent graphite sheets are connected by weak 7r-bonds whereas the 
carbon atoms of a same sheet are connected by cr-bonds. In the 
fullerene the nodeless cr-orbitals are localized at the radius of 
the cage while the single-node 7r-orbitals, on the contrary, are 
oriented perpendicularly to the surface of the fullerene. The ra- 
tio of cr- and 7r-orbitals in C6o should be equal to 3 : 1 due to 
the sp 2 -hybridization of carbon orbitals [22]. 

All 258 electrons of the endohedral compound Ar@C6o are 
described within a framework of the united electronic config- 
uration. The filling scheme is as follows. Double-nodal 3s 2 
shell is added to the configuration of the pristine C6o ©. The 
latter 16 electrons are distributed on the non-closed 10m and 
7h shells so that the number of electrons located on cr- and n- 
orbitals is conserved. As a result the electronic configuration of 
the endohedral compound is written as 

ls 2 2s 2 2p 6 3s 2 3p 6 3d 10 4f 14 5g 18 6h 22 7i 26 8k 30 91 34 1 0m 26 

4d 10 5f I4 6g 18 7h 18 . (2) 

To define single-electron energies and wave functions the 
system of self-consistent equations is solved for all electrons 
of the compound. The potential of positive ions, U, is included 
in the equations and combines the potential of the fullerene core 
and the Coulomb potential due to the nucleus Z of the confined 
atom (Z = 18 in case of Ar), U = [/core - Z/r. Electronic 
system is treated in the HF and the LDA approximations. 

The fullerene core of 60 fourfold charged carbon ions C 4+ 
is described via jellium model as a positively charged spherical 
layer of a finite AR thickness. The average radius R and AR 
are the adjustable parameters of the potential. In the current 
calculations R stands for the radius of Cgo which is equal to 
3.54 A and AR is taken equal to 1.5 A Q. 



3. Results 

The results of the self-consistent calculation of electronic 
structure of Ar@C6o are presented below. Figure Q] represents 
the radial potential of the compound calculated within the HF 
and the LDA frameworks, the 3p and 4d wave functions of 
the compound are presented in Figures [2] and [3] The local 
part of the HF and the LDA radial self-consistent potentials 
are mapped in Figure Q] by the triangles and the filled circles, 
respectively. As far as the system of equations is solved with 
different potentials within the HF approximation and the LDA, 
wave functions obtained within these approximations have a re- 
markable difference. One can see that accounting for the exact 
non-local exchange interaction within the HF approximation re- 
veals significant modification of the 3p and 4d wave functions 
(the line with triangles in Figures [2] and On the contrary, 
accounting for the local exchange interaction within the LDA 
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Figure 1: Radial self-consistent potential of Ai'@Cgo calculated within the HF 
approximation (line with triangles) and the LDA (filled-circled line). Results 
of the LDA calculations performed by Madjet et al. tlal are also presented 
(open-circled line). 



doesn't lead to the such modification of the wave functions (the 
filled-circled line in Figures [2] and[3]l. 

InRef. [16] an assumption on the hybridization of the ar- 
gon 3p shell in Ar@C6o was formulated. In the cited paper the 
system of the LDA equations was solved for the superposition 
of atomic state and the fullerene state (see |24] for the details). 
The hybridized 3p state was defined as a superposition of the 
3p shell of a free argon atom and the 3p shell in pristine Ceo- 
Due to the mixing of the atomic and fullerene wave functions 
the additional node appears in the hybridized 3p function (the 
open-circled line in Figure|2]l. The radial potential of the men- 
tioned calculation is presented in Figure Q] by the open-circled 
line. 
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Figure 2: The 3p wave function in free Ar atom (solid line), in pristine Cgo 
(dashed line) and in Ar@Cgo calculated within the HF approximation (line with 
triangles) and the LDA (filled-circled line). Results of the LDA calculations 
performed by Madjet et al. fT^I are also presented (open-circled line). 
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Figure 3: The 4d wave function in pristine C 60 (dashed line) and in Ar@C 60 cal- Figure 4: Redistribution of the electronic density of the 3p shell in Ar@C 6 o- 
culated within the HF approximation (line with triangles) and the LDA (filled- 
circled line). 



As mentioned above, in the present work we introduce the 
unified electronic configuration. In spite of the remarkable re- 
distribution of the electronic density the 3p wave function has 
a single radial node as opposed to the framework used in Ref. 
Il6ll . Let us note, however, that the symmetry of the wave func- 
tion is significantly changed. The 3p electronic density in the 
pristine Ceo has a minimum at the radius of the fullerene, so it 
corresponds to the 7t-orbital (the dashed line in Figure|2]l. In the 
case of Ar@C6o the 3p shell becomes a cr-type orbital with the 
maximum distribution of the electronic density at the radius of 
the fullerene (the line with triangles). The change of the sym- 
metry and the redistribution of electronic density allows one 
to draw a conclusion on the significant hybridization of the 3p 
shell within the HF approximation. This effect disappears if the 
compound is treated within the LDA due to account of the lo- 
cal exchange interaction between the electrons of the confined 
atom and the delocalized electrons of the fullerene. This result 
generally correlates to the result previously obtained within the 
LDA in lfl6ll (see the open-circled line in FigureUJ). 

For the considered electronic configuration the account 
of the non-local exchange interaction within the HF approxi- 
mation leads to the remarkable charge redistribution due to the 
modification of the 3p and 4d shells (see Figures [4] and |5j. As 
it can be seen in Figure [4] the most part of the electronic den- 
sity of the 3p shell of free Ar is located within the sphere of 
1.5 A. But in the endohedral compound the partial charge of 
the 3p 6 shell equal to 5.7 e (where e stands for the elementary 
charge) is pressed out of this sphere of 1.5 A to the vicinity of 
the fullerene's cage. In contrary, the partial charge of the 4d 10 
shell which is equal to 3.2e is transferred from the vicinity of 
the fullerenes cage to the vicinity of the confined atom. 

4. Conclusion 

To conclude, the self-consistent Hartree-Fock description of 
the endohedral fullerene Ar@C6o was carried out for the first 
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Figure 5: Redistribution of the electronic density of the 4d shell in Ar@Cgo- 



time. Electrons of the encaged atom and the delocalized elec- 
trons of Cgo were considered simultaneously by introduction 
of the united electronic configuration. Having considered the 
configuration described above it was shown that the account of 
the non-local exchange interaction within the Hartree-Fock ap- 
proximation leads to the significant modification of the 3p and 
4d shells as opposed to the local exchange interaction within the 
local density approximation. As a result of the modification the 
redistribution of the electronic density of the 3p and 4d shells 
appears and causes the accumulation of the additional positive 
charge in the vicinity of the encaged atom and the additional 
negative charge on the fullerene core. 

To accomplish in full the analysis of the phenomena of hy- 
bridization and shells modification in Ar@C6o one has to in- 
vestigate electronic configurations other than in (0. This work 
is currently in progress and the results will be published else- 
where. 
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